arXiv:1504.03600v2 [hep-ex] llJun2015 


A Combined Limit on the Neutrino Mass from Neutrinoless Double-/? Decay and 

Constraints on Sterile Majorana Neutrinos 

Pawel Guzowski/ Luke Barnes/ Justin Evans/ Georgia Karagiorgi/ Nathan McCabe/ and Stefan Soldner-Rembold^ 

^ The University of Manchester, Manchester MIS 9PL, United Kingdom 
(Dated: June 12, 2015) 

We present a framework to combine data from the latest neutrinoless double-/! decay experiments 
for multiple isotopes and derive a limit on the effective neutrino mass mpp using the experimental 
energy distributions. The combined limits on mgp range between 130 — 310 meV, where the spread is 
due to different model calculations of nuclear matrix elements (NMEs). The statistical consistency (p 
values) between this result and the signal observation claimed by the Heidelberg-Moscow experiment 
is derived. The limits on mgp are also evaluated in a (3 -I- 1) sterile neutrino model, assuming all 
neutrinos are Majorana particles. 
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I. INTRODUCTION 

The observation of neutrinoless double-/ {Ov/dP) de¬ 
cays would demonstrate the Majorana nature of neu¬ 
trinos [l|, representing direct evidence for physics be¬ 
yond the standard model. Neutrinoless double-/? decay 
is a second-order electroweak process where a nucleus 
decays through the emission of two electrons, {A, Z) 

{A, Z + 2 ) + 2 e~, thereby violating lepton number. The 
half-life of the isotope decaying through Oi// 3/3 decay is 
given by 


2 
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' mj 

assuming the exchange of a light Majorana neutrino 
(mass mechanism). Here, rrie is the electron mass. The 
phase space factors and nuclear matrix elements 
vary with isotope. The decay rate is proportional to the 
square of the effective neutrino mass mpi3, which repre¬ 
sents a coherent sum over the masses rrii of the neutrino 
mass eigenstates, 


m/3/3 


mi \Uel\^ + m2 \Ue2\^ e*“ -I- m 3 |f 7 e 3 |^ 


( 2 ) 


weighted by the squares of the corresponding elements 
Uei of the Pontecorvo Maki Nakagawa Sakata (PMNS) 
matrix and the Majorana phases a and /?. 

The seesaw mechanism is the most common model for 
neutrino mass generation Here, heavy right-handed 
neutrinos mix with the left-handed neutrinos and gener¬ 
ate light Majorana masses for the active neutrinos. In 
addition to the exchange of these light Majorana neu¬ 
trinos, there are several other potential models Q that 
could contribute to Oi// 3/3 decays, e.g., left-right symmet¬ 
ric models, i?-parity violating supersymmetry, or mod¬ 
els with extra dimensions. Observing Oi// 3/3 decays with 
different isotopes could help to disentangle these mecha¬ 
nisms 0. 

The search for OvPP decays is pursued in a range of ex¬ 
periments that use different isotopes. Experiments where 


detector material and Oupp isotope are identical are con¬ 
strained in the choice of isotope by the detector technol¬ 
ogy. This approach is employed by collaborations such as 
GERDA which uses a highmurity ™Ge detector, EXO- 
200 and KamLAND-Zen Q, which contain enriched 
^^®Xe, and CUORICINO [i| and CUORE-0 @, bolometers 
made of Te02 crystals containing ^^°Te. Alternatively, 
detector and isotope can be separated, which allows ex¬ 
periments such as nemo-3 m to study a variety of iso¬ 
topes, e.g., ^™Mo and ®^Se. 

It is therefore important to develop techniques to com¬ 
bine the data taken by different experiments with a range 
of isotopes. This allows us to make quantitative com¬ 
parisons between experiments, to study the consistency 
of their results, and to obtain combined limits on m^/3. 
A previous comparison of different experimental results 
using the published limits on r//2 has been performed 
in Ref. [T^. In this analysis, we simultaneously fit the 
experimental energy distributions to obtain a combined 
limit. 

Such a direct combination of the results for multiple 
isotopes requires a specific NME calculation to relate the 
half-lives of the various isotopes before deriving a limit 
on TO/3/3. We chose a set of commonly used NME mod¬ 
els to perform the combinations, but the procedure is 
generally applicable to include any NME model. The 
combination of the experimental energy distributions is 
performed for each particular NME model. In princi¬ 
ple, different models could give a better estimate of the 
NMEs for different isotopes, but this case is not consid¬ 
ered here, since it would require determining systematic 
correlations between the different models. The final lim¬ 
its on TO/3/3 are given as a range covering the full set of 
NME models, following the procedure normally employed 
by individual experiments. We take into account theo¬ 
retical uncertainties on the calculations where they are 
available. 

In this Article, we present a method based on the pub¬ 
lished energy distributions from all recent experiments 
and use it to derive a first combined limit on TO/3/3 based 
on multiple isotopes. We first demonstrate that we can 
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reproduce the published half-life limits of each individ¬ 
ual experiment, before performing the combination using 
different NME calculations. We study the consistency of 
the combined mpp limit with the positive claim of the 
Heidelberg-Moscow experiment [Q. We also interpret 
the combined limit on m/ 3/3 as a constraint on a (3 -I- 1) 
model with three active and one sterile Majorana neutri¬ 
nos. 


II. METHOD 

Limits are calculated using the CLg method [iHill, 
which applies a modified-frequentist approach using a 
Poisson log-likelihood ratio (LLR) test statistic. The 
value of CLs is defined as the p value of the data un¬ 
der the hypothesis that we observe both signal and back¬ 
ground, CTsb, divided by the p value for the background- 
only hypothesis {CL},). Systematic uncertainties are 
marginalized through Gaussian constraints on the pri¬ 
ors, with the best fits of these parameters determined 
by maximizing the likelihood with respect to the data 
in both the signal-plus-background and background-only 
hypotheses. A limit at the 90% Confidence Level (CL) is 
obtained from the signal strength that produces a value 
of CLs = (100 — 90)%. The expected limit corresponds to 
the median of an ensemble of pseudo-experiments gener¬ 
ated from the background-only probability distributions, 
and the ±1 standard deviation ranges are derived from 
the ensemble. When combining multiple experiments, 
the limit setting process uses a sum of the experiments’ 
individual LLRs, assuming no correlations. 

For a single experiment, the signal strength is inversely 
proportional to the half-life Ty 2 - When combining mul¬ 
tiple experiments, the signal strength is calculated for a 
common mpp, which is related to Ty 2 using Eq. [T] and 
for a specific NME calculation. 

III. EXPERIMENTAL INPUTS 

The most stringent published limits on neutrinoless 
double-/? decay are currently provided by the CUORICINO, 
CUORE, exo- 200, KamLAND-zen, GERDA, and NEMO-3 
experiments. We use their most recent published energy 
distributions, together with the statistical and system¬ 
atic uncertainties on signal and backgrounds, and the 
correlations of the systematic uncertainties as quoted by 
the experiments. The input distributions are shown in 
Fig. [1] with their statistical uncertainties. In total, we 
combine 250 data points. The experimental systematic 
uncertainties are assumed to be uncorrelated between ex¬ 
periments, because the data are taken using completely 
different experimental approaches. 

CUORICINO was a bolometric detector comprising 62 
Te 02 crystals, which were operated between 2003 and 
2008, for a total exposure of 19.75 kg-y Q. Using this 
exposure and the average signal efficiency of 82%, we 


obtain a total expected number of 18 signal decays for 
a half-life of 2.8 x 10^^ years. The energy distribution 
for the signal published in Fig. 9 of Ref. Q normalized 
to this decay rate is used as the input distribution for 
the limit setting. We use an uncertainty for the signal 
normalization of 1.5% and an energy scale uncertainty 
of 0.8 keV. A normalization uncertainty of 4.2% is ap¬ 
plied for the constant background, and the ®°Co back¬ 
ground, which is centred at an energy of 2505 keV, has 
a normalization uncertainty of 7.7%. These background 
uncertainties are the exposure-weighted statistical uncer¬ 
tainties taken from Table 4 of Ref. Q . 

The CUORE experiment has also published results from 
a first phase, CUORE-0, which utilizes a single tower with 
52 Te 02 crystals. The results correspond to an expo¬ 
sure of 9.8 kg-y for the isotope ^^°Te. We use the data 
published in Fig. 3 of Ref. @ and the reported signal 
efficiency of (81.3 ± 0.6)%. The signal shape is taken 
from Ref. [13 , parametrized by a Gaussian function with 
width of « 5 keV. The statistical uncertainty on the 
background normalization is 3.45%, whereas the effect 
of other systematic uncertainties is negligible. 

The exo-200 detector contains a liquid xenon time 
projection chamber that has been in operation since 
2011 A total of 200 kg of xenon is used, with ^^®Xe 
enriched to 81%. We use the energy distribution of the 
single-site (SS) decays published by the exo-200 Collab¬ 
oration in Fig. 4(a) of Ref. @ for an exposure of lOO kg-y, 
since the sensitivity of SS topologies to decays is 

significantly higher than for multi-site (MS) topologies. 
We apply a signal normalization uncertainty of 8 . 6 % with 
an energy scale uncertainty of 0.2%. The background 
normalization uncertainty is 10.9%. We also apply an 
overall normalization uncertainty of 9.6%, correlated be¬ 
tween signal and background distributions, which is de¬ 
rived from the ratio of SS to MS topologies. 

The KamLAND-Zen detector is filled with 13 tons of liq¬ 
uid scintillator, which is loaded with enriched ^^®Xe Q. 
We use Fig. 1(a) of Ref. Q, published after a total expo¬ 
sure 89.5 kg-y with a background normalization uncer¬ 
tainty of 11.21% and a signal normalization uncertainty 
of 3.99%. 

The CERDA Collaboration uses high-purity germanium 
calorimeters enriched in ^®Ge [3. We use the dataset 
obtained between 2011 and 2013 Q, with an exposure of 
21.6 kg-w The energy distribution is taken from Fig. 1 
of Ref. 0 after pulse shape discrimination. We apply a 
9% uncertainty on the signal normalization and a 20% 
uncertainty on the constant background normalization. 

The NEMO-3 Collaboration took data from 2003 to 
2011 with a tracker-calorimeter detector measuring seven 
Ou/3/? candidate isotopes m- Since the ^°°Mo measure¬ 
ments are the most sensitive due to the larger isotope 
mass, we focus on this isotope. The input distribution is 
the data shown in Fig. 2 of Ref. 0, which corresponds 
to an exposure of 34.7 kg-y. The signal normalization 
uncertainty is 7%, and the normalization uncertainty are 
0.7% for the two-neutrino background and 10% for the 
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FIG. 1: Energy distributions for signal and background used as input for the combination with their statistical uncertainties: 
(a) CUORICINO, (b) CUORE-0, (c) EXO-200, (d) KamLAND-zen, (e) GERDA, and (f) nemo- 3 data. The signal distributions are 
normalized to represent the 90% CL observed limits. 


other backgrounds. 

We expect small differences between our method and 
the published results, since we use a different limit set¬ 
ting procedure and only a single input distribution for 
each experiment. As a first step designed to validate our 
method, we determine the lower limits on the half-lives 
at the 90% CL and compare them to the values published 
by the collaborations (see Table |T|. The results for ^®Ge 
and ^°°Mo are in very good agreement with the pub¬ 
lished values. The half-life limit for the exo-200 data 


are 15% — 20% higher with our limit setting method, 
whereas the result for KamLAND-Zen is about 10% lower. 
We observe the same effect for a previously published 
exo-200 result [13, which has also been used by the 
KamLAND-zen Collaboration to set a combined limit for 
^^®Xe (3- We combine the same two ^^®Xe data sets and 
obtain a limit of = 34 x 10^"^ y, in perfect agreement 
with the published KamLAND-zen combination. 

We also obtain perfect agreement for the CUORICINO 
result using ^^°Te, while the limit on obtained by the 
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Limit on 

-^1/2 

( 10 ^" y) 


Experiment 

Publ. 

Obs. 

Exp 

. ±1(7 range 

1-CLb 

CUORICINO 

2.8 

2.8 

2.9 

2.0 - 4.2 

0.474 

CUORE-0 

2.7 

3.0 

3.0 

2.1 - 4.3 

0.520 

Combined 

4.0 [9] 

4.4 

4.3 

2.9 - 6.2 

0.513 

exo-200 

11 m 

13 

21 

14 - 30 

0.131 

KamLAND-zen 

19 [7] 

17 

11 

7-15 

0.918 

Combined 

— 

21 

24 

16 - 34 

0.360 

GERDA 

21 [a 

20 

21 

14 - 29 

0.450 

nemo-3 

1.1 fll] 

1.1 

0.9 

0.6 - 1.4 

0.634 


TABLE I: The published limits on for each experiment 
are compared to the calculated observed and expected limits. 
The ±1(T range around the expected limit and the 1 — CL^ 
value of the data are also shown. 

CUORE-0 Collaboration is 10% lower than the limit we 
derive. Since CUORICINO and CUORE-0 are two phases of 
the same experiment, we only use the combined CUORE 
limit to perform our remaining analysis. 

The observed limits are also compared to the median 
expected limits and the corresponding ±lcr range around 
the median expected limits. All observed limits lie within 
the ±1(7 band, except for the KamLAND-Zen result, where 
the observed limit of = 17 x 10^"^ y is about 1.5 (t 
better than the expected limit, corresponding to 1 —CL^ 
of 0.918. This is caused by the dehcit of data compared 
to the background expectation in the signal region. 


IV. NUCLEAR MATRIX ELEMENTS 

We derive limits on the effective neutrino mass for dif¬ 
ferent NME calculations, which are summarized in Ta¬ 
ble EH The five different models used in the combina¬ 
tion are the Generating Coordinate Method (GCM) 
the Interacting Boson Model (IBM-2) [I^, the Inter¬ 
acting Shell Model (NSM) [2^, and three models using 
the Quasiparticle Random-Phase Approximation: (i) the 
QRPA model using four different sub-calculations based 
on the Argonne V18 (A) or charge-dependent Bonn (B) 
nucleon-nucleon potentials with an older (o/d) or newer 
(new) version of the parametrization of the particle- 
particle interactio ns l^ . (ii) the proton-neutron QRPA 
model (pnQRPA) , and (hi) a QRPA model that uses 
its renormalized version for evaluating uncertainties, la¬ 
belled (R)QRPA m. 

The model calculations predict values for the NMEs 
which can differ by up to a factor of ~ 2 for a particular 
isotope. The NSM model predicts the smallest NMEs for 
all isotopes, apart from ^°°Mo which is not evaluated in 
this calculation, whereas the models differ in their pre¬ 
dictions of the isotopes with the largest NMEs. 

The (R)QRPA and the IBM-2 calculations also provide 


uncertainties on the NMEs. The IBM-2 model uncertain¬ 
ties on the NMEs are 16%, which we have assumed to be 
fully correlated across isotopes. The (R)QRPA model in¬ 
cludes a correlation matrix between isotopes, which we 
take into account in the limit setting. The uncertainties 
are evaluated in Ref. by using (i) two values of the 
weak axial-vector coupling parameter, = 1 and 1.25, 

(ii) two different approaches to short-range correlations, 

(iii) two many-body models, QPRA and RQPRA, and 

(iv) three different sets of single-particle states. The un¬ 
certainties also include the experimental uncertainties on 
the particle-particle coupling constant gpp extracted from 
2 jz/ 3,5 data. The resulting 24 combinations are used to ex¬ 
tract correlation coefficients from the error ellipses which 
cover the full range of possible outcomes. Eor our exper¬ 
imental combination, we reinterpret these model uncer¬ 
tainties as Gaussian uncertainties by reducing them by 
a factor of 0.68. We quote the reduced uncertainties in 

Table EH 

The phas e sp ace factors are taken from a recent 
calculation [^. We set the weak axial-vector coupling 
parameter to gA = 1.27. The uncertainties on the phase 
space factors originate from the uncertainties on the Q 
values and the nuclear radius. The uncertainty on the 
nuclear radius dominates and leads to an uncertainty on 
the phase space factor of w 7%. 


V. RESULTS 

Only in the case of ^^®Xe, two independent experi¬ 
ments, KamLAND-Zen and exo- 200, measured the same 
isotope. Their data can be combined directly without 
using NMEs to derive a limit on Ty 2 - combination 
yields a limit of Tjy '2 > 2.1 x 10^® y, which is dominated 
by the exo- 200 measurement due to the higher exposure 
compared to KamLAND-Zen. This observed value is less 
stringent than the previously published combined value 
of T^J 2 ^ ^ y 0 based on a smaller exo- 200 data 

set [^, which displays a downward fluctuation of the 
data relative to the background expectation. The com¬ 
bined sensitivity, given by the median expected limit, has 
improved from = 1.6 x 10^® y to = 2.4 x 10^® y 
due to the increased exposure. 

To combine results for multiple isotopes, we include 
the effects of NMEs. We show an example in Fig. [2] for 
the GCM calculation. The results for mpp are shown 
for each individual experiment, ordered by their sensi¬ 
tivities using this model. The combined observed limit is 
TO/ 3/3 < 130 meV for an expected limit of TO/ 3/3 < 130 meV, 
which improves upon the best individual limit obtained 
by KamLAND-Zen of mpp < 150 meV and the best indi¬ 
vidual expected limit by exo-200 of mpp < 140 meV. 
All limits are given at the 90% CL. Since the mass limit 
depends on the fourth root of the exposure, these im¬ 
provements corresponds to an increase in exposure by a 
factor of 1.4-1. 6 . 
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Isotope 


Phase Space 
Factor G°'' 

(lo-^v) m 


Nuclear Matrix Element Models 

GCM IBM-2 NSM QRPA pnQRPA [22] (R)QRPA [2^ 

[la [1^ [20] A-old A-new B-old B-new NME Rel. Unc. Correlation Matrix 



0.615 

4.60 

4.68 

2.30 

5.812 

5.157 

6.228 

5.571 

5.26 

4.315 

0.191 

1 


i°°Mo 

4.142 

5.08 

4.22 

— 

5.696 

5.402 

6.148 

5.850 

3.90 

3.184 

0.254 

0.973 

1 

isoTe 

3.699 

5.13 

3.70 

2.12 

4.306 

3.888 

4.810 

4.373 

4.00 

3.148 

0.247 

0.899 

0.862 1 

i36Xe 

3.793 

4.20 

3.05 

1.76 

2.437 

2.177 

2.735 

2.460 

2.91 

1.795 

0.293 

0.805 

0.747 0.916 1 


TABLE II: Phase space factors for qa = 1.27 for the four isotopes, values of the nuclear matrix elements for the GCM, IBM-2, 
NSM, QRPA, pnQRPA, and for the (R)QRPA NME calculation, together with the relative uncertainties on the (R)QRPA 
NMEs and their correlation matrix. The relative uncertainties quoted for the IBM-2 calculation are 0.16 for each isotope. 



m (meVj 


FIG. 2: Individual experiment effective mass limits, and the 
combined limit, using the GCM model. 

Using the GCM model, the CUORE and GERDA data 
have sensitivities of mpp = 250 meV and 310 meV, re¬ 
spectively, while the two experiments using ^^®Xe have 
the best sensitivities in the range mpp = 140-190 meV. 
Figure 13] shows the observed and expected limits for each 



FIG. 3: Observed and expected limit onm^p for the different 
experiments and NME models without NME model uncer¬ 
tainties. 

NME model and experiment separately. For this compar¬ 
ison, NME uncertainties are not taken into account, since 


only some of the calculations provide them. The ordering 
of the experimental sensitivities changes with the NME 
model used to derive the mpp limit, further emphasising 
the need for using a range of isotopes to search for Qvfij3 
decays. 

The coefficients of the correlation matrix provided by 
the (R)QRPA calculations for the four isotopes consid¬ 
ered here range from 0.747 to 0.973 (see Table P. To 
study the effect of these correlations, we derive the limits 
on mpp assuming no NME uncertainties, the full cor¬ 
relation matrix, and a single correlation coefficient for 
all isotope combinations ranging from —0.3 to 1. The 



FIG. 4: Limit on mpp as a function of a single correlation coef¬ 
ficient p in the (R)QRPA model. The lines at mpp = 250 meV 
and 290 meV represent the limits without uncertainties and 
with the full correlation matrix of the (R)QRPA model, re¬ 
spectively. The red line shows the range of correlation coeffi¬ 
cients. 

results are shown in Fig. |4l Using no NME uncertain¬ 
ties yields a limit of mpp < 250 meV, which increases 
to mpp < 290 meV once the model uncertainties and 
the full correlation matrix are taken into account, and to 
mpp < 300 meV assuming 100% correlation between the 
NME uncertainties for the different isotopes. The vari¬ 
ation in limits of Ri 40 — 50 meV is comparable to the 
range observed for the different model assumptions with 
the QPRA calculation, which are expected to represent 
similar model uncertainties. 

In Table Ell we show the combined limits on mpp for 
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NME 

.™obs 

^/3I3 

Improvement 

p value 


(meV) ( 

meV) 

Limit 

Sensitivity 

(HM) 

GGM 

130 

120 

14% 

(K) 

10% (E) 

0.001 

IBM-2: 







no unc. 

170 

170 

16% 

(K) 

12% (E) 

0.010 

with unc. 

190 

180 

16% 

(K) 

13% (E) 

0.021 

NSM 

310 

290 

14% 

(K) 

10% (E) 

0.003 

QRPA: 







A-new 

200 

200 

23% 

(G) 

25% (E) 

0.095 

A-old 

180 

180 

26% 

(G) 

25% (E) 

0.100 

B-new 

180 

180 

28% 

(K) 

24% (E) 

0.073 

B-old 

170 

160 

28% 

(K) 

23% (E) 

0.077 

pnQRPA 

170 

170 

19% 

(K) 

16% (E) 

0.029 

(R)QRPA: 







no unc. 

250 

240 

25% 

(G) 

25% (E) 

0.109 

with unc. 

290 

290 

23% 

(G) 

21% (G) 

0.311 


TABLE III: For each NME calculation, the combined ob¬ 
served and expected limits on mpp, the improvement in the 
limit and the sensitivity relative to the best individual experi¬ 
ments for that NME model. The best experiments are gerda 
(G), exo-200 (E), or KamLAND-zen (K). The p value of the 
limit with respect to the Heidelberg-Moscow (HM) positive 
claim [I^ are also shown. All mpp limits are given to two 
significant digits. 



FIG. 5: Mass limits (observed and expected) for each NME 
calculation. For the IBM-2 and (R)QRPA models, the limits 
are given both with and without the uncertainty on the NME 
calculation. The effective mass for the positive Heidelberg- 
Moscow (HM) claim [l^ is also shown. 


each NME model. We provide the observed and expected 
limits for each combination and the improvements in ob¬ 
served limit and sensitivity over the result obtained by 
the best individual experiment. As shown in Fig. [3l the 
most sensitive experiments change for the different mod¬ 
els. 

The combined results are compared to the result ob¬ 
tained by the Heidelberg-Moscow (HM) experiment in 
Fig.El where the error bands represent the uncertainties 
of the HM half-life measurement of T^j 2 ~ ( 2 - 231 ' 1 q' 3 i) x 
10^® y for ^®Ge [l^. The p values, testing the consistency 


between the result obtained by the HM experiment us¬ 
ing ^®Ge and the combined upper limits, are also given 
in Table IHII They are calculated taking into account the 
experimental uncertainties. 

The p values differ significantly between different NME 
models, and any conclusion on the level of consistency be¬ 
tween the combination of the latest experimental results 
and the result obtained with the the T^j 2 reported by the 
HM group depends strongly on the chosen NME model. 
The p values range from 0.001 (about 3 standard devi¬ 
ations) for the GGM model to 0.311 for the (R)QRPA 
models. The p value for the gerda result and this HM 
measurement is 0.63, independent of NME model since 
the experiments both use the isotope ^®Ge. 



FIG. 6: Gombined limits on the effective mass mpp as a func¬ 
tion of the lightest neutrino mass mo. The width of the hor¬ 
izontal band represents the range of limits obtained with the 
different NME models. The neutrino oscillation parameters 
are taken from Ref. [^, using the best fit values. The inner 
band for the NH and IH fits shows the full range associated 
with the unknown Majorana phases and the outer band the 
additional effect of the 3-standard deviation experimental un¬ 
certainties on the oscillation parameters. The limit on mo 
derived from the Planck data is also shown. 

The best published limit on mpp by KamLAND-Zen Q 
is obtained by combining their data with the earlier 
exo-200 result 0. Using a similar set of NMEs, the 
mass range they obtain is mpp = 120-250 meV, whereas 
our observed combined limit in the range mpp = 130- 
310 meV is slightly higher because we use the latest exo- 
200 data set. 

In Fig. [6l we overlay the range of observed combined 
limits on the allowed effective mass mpp as a function 
of the lightest neutrino mass toq, with the bands rep¬ 
resenting current measurement uncertainties of the neu¬ 
trino oscillation parameters assuming either the normal 
hierarchy (NH) with mi < < mp or the inverted hi¬ 

erarchy (IH) with mp < mi < m 2 for the ordering of the 
neutrino mass eigenstates in the standard three-neutrino 
framework. We also show a limit on mg, which is derived 
from the limit on the sum of the neutrino masses ob¬ 
tained by the Planck Gollaboration [ 2 ^ . They set a limit 









































7 


of rrii < 492 meV at the 95% CL based on the tem¬ 
perature and the polarization of the Cosmic Microwave 
Background. 

We also interpret the combined limit on as a con¬ 
straint on a single sterile Majorana neutrino that mixes 
with the three active Majorana states. In this model, the 
expression for the effective neutrino mass defined in 
Eq. [2] is modified to take into account the mixing with 
the fourth neutrino state: 


771/3/3 = 


mi|17el|"+77l2|17e2|"e*“ + 
+m3 pesf + rn4 |C/e4|^ e*'*' 


(3) 


It includes an additional mass term 7774 , a Majorana 
phase 7, and the element Uei of the now extended, 
4x4 mixing matrix, which is related to the sterile neu¬ 
trino mixing angle ^14 through |C/e 4 |^ = sin^ @14 in the 
parametrization used here. In this model, the unitarity 
constraint applies only for the 4x4 mixing matrix. How¬ 
ever, the central values of the elements Uei{i = 1, 2 , 3) do 
not change significantly when the 3x3 unitarity con¬ 
straint is removed from the global fit [l^. The unitar¬ 
ity constraint on the 4x4 mixing matrix will restrict 
|C/e 4 p = sin^ 014 to be < 0 . 1 . 

Our ability to probe a certain range of 777.4 |C^e 4 |^ de¬ 
pends on 7770 and the Majorana phases (see, e.g., (^[^ 1. 
We take the approach of translating the limits on 777 / 3/3 
into corresponding limits on the combination 7774 |C/e 4 |^ 
under two extreme situations for a given 7779 : 

• The NME model (NSM) that predicts the highest 
(least stringent) limit on mpp, together with the 
phases a = j5 = 0 , and 7 = tt that provide the 
smallest contribution to 777 / 3 / 3 . 

• The NME model (GCM) that predicts the lowest 
(most stringent) limit on 777 / 3 / 3 , together with the 
phases a = /3 = 7 = 0 that provide the largest 
contribution to mpp. 

This first parameter set allows for large contributions 
from 7774 |17e4|^ and thus provides the least constraining 
limit on 7774 |17e4|^, whereas the second set allows for only 
small contributions from 7774 |17e4|^ and thus yields the 
most constraining limit on 7774 |C/e 4 |^. 

In Figs. Hand|S1 we translate the measured combined 
limit on 777^/3 for the two situations above into a limit in 
the ( 7774 , |t/e 4 |^) plane for a set of different assumptions 
on 7779 . To be consistent with cosmological bounds, we 
vary mg between 0 and 150 meV. For any of the 7779 
values considered, all other NME model and Majorana 
phase combinations will produce limit curves that lie in 
between these extremes. The values of and CAi 

{i,j = 1,2,3) are taken from a recent global fit [^ . 
Including the three-standard deviation uncertainties on 
these fit values has negligible effect on the results. 

The central value of sin^ 20i4 = 0.09 and = 

1.78 eV^ from a global fit of the oscillation parame¬ 
ters in a (3 -|- 1) sterile neutrino model including reactor. 



FIG. 7: Limits on mpp for different NME models translated 
into a constraint on the sterile neutrino in the (7774, | 17 e 4 |^) 
plane in a (3 -I- 1) model. Different values of the Majorana 
phases a, P, 7 and of the lightest active neutrino mass mi 
in the NH are shown. The complete region is excluded for 
nil — 150 meV with the GCM model and vanishing Majorana 
phases. 



FIG. 8 : Limits on for different NME models translated 
into a constraint on the sterile neutrino in the (m 4 , \Uei\'^) 
plane in a (3 -I- 1) sterile neutrino model. Different values of 
the Majorana phases a, /3 ,7 and of the lightest active neu¬ 
trino mass m 3 in the IH are shown. The complete region is 
excluded for m 3 = 150 meV with the GCM model and van¬ 
ishing Majorana phases. 


gallium, solar neutrino, and lsnd/karmen v,, disappear¬ 
ance data, all of which are directly sensitive to 7774 and 
|C/e 4 |, is also shown in Figs. [7] and [5] The limit from the 
077/3/3 combined fit is currently close to the (3-1-1) global 
best fit value only in the more constraining scenarios cal¬ 
culated using the GCM model and with 7779 > 100 meV 
and favorable Majorana phases. Our results are consis¬ 
tent with previous analyses (see, e.g., [H, [3l|) that in¬ 
vestigate 777^/3 regions allowed by global fits to oscilla¬ 
tion experiments under a (3-1-1) hypothesis. Assuming 
the sterile neutrino is a Majorana particle, neutrinoless 










double-/? decay can therefore independently constrain the 
existence of a sterile neutrino under those extreme as¬ 
sumptions. For a more detailed discussion on constraints 
to (3-1-1) sterile neutrino models from neutrinoless double 
/? decay and projected implications for next-generation 
experiments, see Ref. (^ . 

VI. SUMMARY 

We have performed the first combination of the lat¬ 
est data sets from experiments searching for neutrinoless 
double-/? decay using multiple isotopes. Using the CLg 
method, we set a limit on the effective neutrino mass 
mpp in the range 130 — 310 meV, depending on NME 
model. Combining the data from multiple isotopes and 
experiments significantly increases the sensitivity over us¬ 
ing the single-best experiment only, corresponding to an 


increase by a factor of Ri 1.5 — 2.4 in exposure. We 
compare these limits with the claimed observation of the 
Heidelberg-Moscow experiment and obtain p values that 
differ significantly depending on the NME calculations 
chosen, ranging from 0.001 for the GCM model to 0.31 
for the (R)QRPA model. Using the uncertainties and the 
full correlation matrix provided by the (R)QRPA model 
changes the limit on mpp by 40 — 50 meV compared to 
using no NME uncertainties. We also translate the com¬ 
bined limit on mpp into a constraint on a light sterile 
Majorana neutrino in a (3 -I- 1) model. This translated 
limit is NME and uiq and Majorana phase dependent. 
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